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I. INTRODUCTION 

The observations by the WMAP satellite HSU, HQ 
have brought the global cosmological data set up to a 
quality where, for the first time, it is possible to obtain 
precision constraints on cosmological models. That the 
data provides no indication of any significant departure 
from gaussianity, adiabaticity, or scale-invariance, and 
furthermore reveals two coherent peaks in the spectrum 
of cosmic microwave background (CMB) anisotropics, 
lends powerful support to the idea of inflation in the early 
Universe as a source for the observed perturbations. This 
opens the prospect of constraining and excluding regions 
of inflation model parameter space. 

The qualitative breakthrough of the WMAP data is 
important. For the first time we have available a CMB 
spectrum that spans from cosmic variance limited mea- 
surements on large angular scales across to the measure- 
ment without any overall calibration error of the peaks on 
small angular scales. Measurements of the polarization 
of the CMB |(| provide some insight into the epoch of 
reionization which in turn helps to constrain inflationary 
models by limiting the effects of parameter degeneracies. 

In this paper we analyze slow-roll inflation models, fol- 
lowing the strategy outlined in Leach et al. 0, and pro- 
visionally applied to pre- WMAP CMB data by Leach 
and Liddle y| as a demonstration of methodology. With 
WMAP it is possible to make the first serious applica- 
tion. As compared to our earlier work, we make sev- 
eral improvements. We use the now-ubiquitous Markov 
Chain Monte Carlo method 0, to obtain the like- 
lihood function over parameter space, we include both 
short-scale CMB data and the galaxy power spectrum 
data from 2dF (but not any lyman-alpha data, whose in- 
clusion has proven controversial |ll|). and we study the 
effect of varying one further slow-roll parameter. Our 
approach differs in several respects from the other pa- 
pers that have already appeared discussing inflation post- 

wmap [unm, 

and we contrast our work with theirs 

in the conclusions. 



II. METHODOLOGY 

We follow the methodology described in Refs. 0, @, 
to which the reader is referred for further details. In 
terms of the horizon- flow parameters ei, ti, etc., the in- 
flationary scalar and tensor power spectra can be well 
represented as power-laws with amplitude and spectral 
indices given by 0, 0] 

A s = H 2 (l-2(C+l)ei-Ce 2 ) (1) 
7reimp L 

n s -l = -2ei-e 2 (2) 

n T = -2ei, (3) 

where C ~ —0.73. The relative amplitude of tensor and 
scalar perturbations is given by 

R = 16ei . (4) 

Later we will also consider weak running (scale- 
dependence) of the scalar spectral index. Although the 
full inflationary predictions are somewhat more detailed, 
Eqs. capture the essence of the inflationary power 

spectra. In the following analysis we start by using 
the first-order power-law shape predictions for inflation 
which include the Stewart-Lyth correction to the am- 
plitude, but ignore any term 0(e 2 ). Later we use the 
full second-order predictions where the terms 0(e 2 ) are 
included in the fit. 

The data that we use in this paper comes from 
VSA O, CBI [lp, ACBARO, WMAP @ and the 
2dF galaxy redshift survey We compute the mi- 

crowave anisotropies using the CAMB code [2(1 cou- 
pled to our own slow- roll inflation module 7], and 
use the package CosmoMC lhC|. modified to include 
the WMAP likelihood code |2l|, in order to com- 
pute the likelihood over parameter space. We gener- 
ate a Markov chain of 60,000 elements. We assume 
a flat ACDM universe and adopt the parameter ba- 
sis {w B , w D , H , 10 W A s (k*), z ro, ei(fc*), £2(fc*)} where ujb 
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FIG. 1: ID posterior constraints for the basic cosmological 
parameters assuming slow-roll inflation. 



and wd are the baryon and dark matter physical den- 
sities, Hq is the Hubble constant, As is the ampli- 
tude of scalar perturbations, z te is the redshift of reion- 
ization (which is assumed to be instantaneous), and 
fc* = O.OlMpc -1 . The only prior that has any effect 
on the constraints is insisting z ro > 4. We are also 
interested in the constraints on the derived parameters 
{Ase~ 2r , ng — 1, -Rio} where r is the optical depth to the 
last scattering surface and Rio 
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III. THE CONSTRAINTS 
A. Constraints on slow-roll inflation 

We begin by considering the constraints on slow-roll 
inflation models, initially only including the horizon- 
flow parameters e\ and ti- For orientation and 
comparison with other work, in Fig. ^ we display 
the constraints on the basic cosmological parameters 
{wb, wd, H , f0 10 As(fc*), Zrc}- Our results are in good 
agreement with other authors, unsurprisingly as we use 
many of the same codes and a similar data compilation. 

Fig. |21 shows the likelihood distribution in the plane 
of the horizon- flow parameters e± and ti- We can see 
that the constraint on €2 is highly correlated with ei, 
since both parameters contribute to the spectral index. 
Moreover the data introduce a further degeneracy, the 
tensor degeneracy. This occurs where models with a ten- 
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FIG. 2: 2D posterior constraints in the £1-62 plane, 
contours are the lcr and 2a bounds. 



The 



sor component, and hence more power on large scales, 
require more power to short scales, and hence a bluer 
spectrum. This is clear from Fig. [21 in which we plot the 
same constraints in terms of the derived parameters ns 
and Riq. 

Because of the strong degeneracy between t\ and 62, 
we read off the 2a upper limit on e\ from Fig. [2] without 
marginalizing out £2, finding 



61 < 0.032, 



(5) 



which gives the best measure of the relative (primordial) 
contribution of tensors, via Eq. Q . This constraint is in 
agreement with Refs. UE]- The direct contribution of 
the tensor spectrum to the Cg spectrum is also of some 
interest and, similarly, we obtain the 2a upper limit 



i?io < 0.32, 



(G) 



and using Eq. we obtain an upper limit on the energy 
scale of inflation 



< 1.4 x 10 -5 . 

m P1 



(7) 



The constraints on the horizon-flow parameters are best 
summarized in Fig. [3 However, we can dehne a new pa- 
rameter along the tensor degeneracy direction and obtain 
the constraint 



£2 - 

10 10 A s e- 2T 



5.0ei 
-82ei 



0.036 ±0.025, 
19.3 ±0.7, 



(8) 
(9) 
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FIG. 3: 2D posterior constraints in the (n B — l)-Rio plane, 
again at la and 2a. Models with a tensor spectrum on large 
scales require a bluer scalar spectrum in order to increase 
CMB power to short scales. 



and these constraints are displayed in Fig. 0] To obtain 
constraints on the shape of the inflaton potential we use 
the slow-roll approximation 
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and rewriting the constraints Eq. ©-© we find 

V 



< 0.23 x 10" 
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(10) 

(11) 
(12) 

(13) 



pi 



'"pi 
4tt 



16tt ^ V 

yl \ 2 y„ 

2.25 | — - — 
V V 



< 0.032, (14) 
= 0.036 ±0.025. (15) 



From these constraints and from Fig. [3 we see that there 
is no evidence of deviations from the extreme slow-roll 
limit Ci — 0, which lies comfortably within the la con- 
tour. The best one can say is that the tensors and tilt 
can be used to marginally improve the fit to the data. 

In comparison with the situation before WMAP 
(e.g. Ref. (3), the principal change is a considerable tight- 
ening of the uncertainties, with the general trend of the 
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FIG. 4: ID posterior constraints on inflationary parameters. 
The solid line corresponds to a power-law fit to the data using 
ei and £2. The dashed line corresponds to a fit where weak 
running of the spectral index is included in the fit via the 
slow- roll parameter £3, which is unconstrained by the data. 



shrinking being towards the scale-invariant case. This is 
sufficient to exclude a significant chunk of slow-roll infla- 
tion parameter space. Of the models remaining, there is 
a mild preference for a red spectral index (n < 1) but not 
with any significance. Because of the tight correlation be- 
tween t\ and €2 there exists a class of models with blue 
scalar spectra which show a mild preference for a tensor 
component on large scales, as is visible in Fig [3] This is 
an area of parameter space which has not been signifi- 
cantly populated with models, as models with blue spec- 
tra tend to have negligible tensors. These models have 
£2 < meaning that the fractional kinetic energy of the 
inflaton is decreasing, which could correspond to models 
leaving a kinetic energy dominated epoch. These models 
are also somewhat 'protected' by the tensor degeneracy 
direction and could prove quite resistant to observational 
pressure for some time to come. 

We now ask how well motivated it was to stop at 62, 
by including one further horizon- flow parameter £3. This 
allows us to include a running (scale-dependence) of the 
scalar spectral index, given by 



din k 



-2eie 2 - e 2 e 3 



(16) 



and to take the expressions for the spectral indices them- 
selves to second-order [l^. The effect of including this 
extra parameter is shown in Fig. 0] where it has a fairly 



4 



modest effect on the fikelihood distributions for other 
inflationary parameters. However e 3 itself is poorly con- 
strained, and is readily consistent with zero. We conclude 
that there is clearly no motivation to include this extra 
parameter, with the improved goodness-of-fit being insuf- 
ficient to warrant its inclusion. We return to this issue 
in Section HVBl 



B. Constraints on power- law inflation 

If we have a specific class of inflation models in mind 
then we can go beyond the constraints of Fig. [3 and 
in this section we examine the constraints on power-law 
inflation. This remains an interesting model because the 
potential, once normalized to the observed perturbation 
amplitude, is described by a single parameter, and so we 
can expect tighter constraints than in the case of general 
slow-roll models. Power-law inflation is expansion 
given by 



a oc t p , p > 1, 

ex = -, e l = 0, i > 2 . 
P 

Equivalently, in terms of conformal time we have 

a oc \i]\ q , — oo < q < — 1 , 
2 1 

<7 



l + 3w 



1-ei 



(17) 
(18) 

(19) 
(20) 



We obtain the constraint on power-law inflation by read- 
ing off the 2cr bound on e\ at the intersection with the e\ 
axis and find: 



< ei < 0.019, 
p > 53, 
-1.019 <g<-l. 



(21) 
(22) 
(23) 



The constraint on e\ is tighter than for slow-roll models 
because power-law inflation requires a red scalar power 
spectrum, and so there is no possibility of taking advan- 
tage of the tensor degeneracy. If in the near future the 
HZ spectrum is ruled out in favor of red-tilted spectrum 
then we can begin to place an upper limit on the index 
p, and needless to say, power-law inflation can be ruled 
out altogether if a blue-tilted spectrum is favored or if 
the tensor spectrum has the wrong relationship with the 
scalar spectrum. 

An easy way to overinterpret the data would be to 
point out that our best-fit models lie close to the region 
well described by power-law inflation (ie. the e% axis), 
with inflationary parameters t\ ~ 0.01 and 62 ~ 0.00. 
However this is likely to be a result of the mild preference 
for red scalar spectra combined with a small slide along 
the (flattish) tensor degeneracy direction once tensors are 
included in the fit. Nonetheless, it is intriguing that such 
a simple inflation model should provide such an excellent 
fit to the data. 



C. Constraints on monomial inflation 

In this section we examine the constraints on monomial 
inflation potentials of the form 



\m P iJ 



(24) 



Once normalized to observations these models are defined 
by two parameters, the index a and 0*, the value of <f> 
when the scale fc* crossed the horizon during inflation. 
However, there is an extra ingredient: the energy scale 
of inflation is fixed once we specify As and ei, which is 
typically around 10 16 GeV. This in turn specifies the max- 
imum number of e-folds of slow- roll inflation after horizon 
scale crossing, iVj^ x , after which inflation must end (re- 
gardless of the mechanism that ends inflation) giv ing way 
to reheating and the standard expansion history |23L |24| . 
In addition, monomial inflation provides a mechanism to 
end inflation via the violation of slow-roll, and so we can 
calculate the functions H(N), ei(N), and 62(N) where N 
is the number of e-folds from the end of inflation. This 
allows us to map a constraint on iVhor to a constraint on 
the slow-roll parameters. 

Substitution of the potential of Eq. lf^|l into Eq. l(TT|) 
and l|12f> gives 



ei 



(2 



167T \TOpi 



a / cj> \ _ 4 
47r \mp\ J a 



(25) 
(26) 



We can calculate the number of e-folds of inflation from 
the definition 



dN 

^ N 



dei 

1 {N) £ 1 £ 2 ' 



(27) 
(28) 



where ei(N) is the initial value of t\. Using 62 = 4ei/a, 
which is a good approximation for monomial inflation, 
we have 



N=™ 
4 



1 

1 

ei 



(29) 



and in the limit N 3> a/4 the horizon flow parameters 
are given by 



a 

AN' 
1 

Eq. H31|) gives us the useful constraint 

1 

£2 > 



(30) 
(31) 

(32) 
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FIG. 5: 2D posterior constraints in the £1-62 plane, for the 
region £2 > 0. The contours are the la, 2a and 3a bounds. 
The hatched region £2 < 1/60 is inaccessible to monomial 
inflation models. The thick lines indicate the available pa- 
rameter space for two monomial inflation models: the Xcj) 4 
model is under strong pressure from observations. 



inflation index is constrained from above to be 

a < 4.3, 3cr 
a < 3.9, 2a 

a < 2.8, lcr. (33) 

Thus, while it is still too early to definitively rule out 
the \<f) A inflation model, it is clear from the pattern of 
the above constraints that this particular model is un- 
der very strong pressure from observations. The m 2 <f> 
is valid as long as 32 < iVhor < 60 (note that monomial 
inflation models with iVh or < 55 must have a prolonged 
reheating epoch at the end of inflation or non-standard 
post-inflationary evolution [3). It is clear from Fig. [5] 
that the same kind of pressure will build on the m 2 (j) 2 
model if the spectral index becomes more tightly con- 
strained, since lines of constant spectral index correspond 
to ei = — (f2 + — l)/2, which, modulo the tensor de- 
generacy, runs parallel to the constraint contours. 

It is natural to ask what can we say in a model- 
independent manner about slow-roll inflation models 
while still including the bound on -/Vhor- The number of e- 
folds of inflation, given by Eq. depends only on the 
function 62(61), and this gives us a straightforward way 
to study the phenomenology of inflation models. Unfor- 
tunately, we only have information about the function 
£2(61) across the 8 or so e-folds across observable scales, 
and so the rest of this function can be at best extrapo- 
lated from observable scales. The initial slope is given by 



where 2V* is the number of e-foldings from the end of 
inflation that the scale fc* crossed the horizon (in the 
language of Ref. 0, iV» ~ iV h or - 4). It should be em- 
phasized that the specific form of this constraint applies 
only for the monomial potentials. The maximum number 
of e-folds of slow-roll inflation corresponding to our scale 
fc* = 0.01 Mpc -1 is approximately 60 |23l l24l | and so we 
have the constraint £2 > 0.017. Using this additional con- 
straint we can rule out certain monomial inflation models 
independent of the actual number of e-folds of slow-roll 
inflation after horizon scale crossing. 

In Fig. EI which is a zoom of part of Fig. [3 we restrict 
our attention to £2 > which corresponds to models 
where the ratio between the scalar field kinetic and total 
energy density is increasing, arguably the most natural 
candidates for ending inflation by the violation of slow- 
roll. The solid lines show the location of the quadratic 
and quartic potential models at different numbers of e- 
foldings up to the maximum. The quartic potential lies 
outside the 2a contour for any allowed value of iV». 

We can constrain the exponent a by reading off the 
value of 4ei x jV™ ax at the point where the observational 
constraints on the slow-roll parameters intersect with the 
constraint on £2 given by Eq. I.'i2[l . We find the monomial 



7T = * ( 34 > 

«£2 £3 

and so the higher slow-roll parameters, via the running 
of the spectral index Eq. i|16|) , would contain useful addi- 
tional information about the model of inflation, assum- 
ing inflation occurs at these high energies. Since we have 
found no constraint on £3 , then it is clear that we can not 
make this type of model independent analysis at present. 
Note also that the bound JV* < 60 can be relaxed some- 
what to + AN if H is significantly reduced in the late 
stages of inflation [24| • In terms of the slow-roll parame- 
ters 

LN ~f * (35) 

For instance, this correction is relevant for the A</> 4 which 
has AN ~ 4, although this is certainly not enough to take 
much pressure off this model. 

IV. DISCUSSION 

A. The Harrison— Zel'dovich spectrum 

We have discounted the inclusion of £3 because the data 
are unable to distinguish it from zero. However, following 
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FIG. 6: ID posterior constraints for the basic cosmological 
parameters. The dotted line corresponds to the Harrison- 
Zel'dovich model and the full line corresponds to inflationary 
models (as shown in Fig. 0. 



the logic of that statement we are forced to conclude 
that neither t\ nor ti have been shown to be inconsistent 
with zero. The minimal assumption concerning the data 
is in fact that it is due to a Harrison-ZePdovich (HZ) 
spectrum, with no indication of inflationary dynamics. 



It is plausible therefore to also regard t\ and €2 as giv- 
ing poor value, in terms of the improvement to the fit 
they give, for their inclusion as extra parameters. It is 
therefore interesting to see how constraints tighten up if 
we assume an HZ spectrum in place of slow-roll infla- 
tion. In Fig. we display once more the constraints on 
the basic cosmological parameters, but now also showing 
the constraints assuming the HZ spectrum and no ten- 
sors. It is clear that, broadly speaking, the two models 
are in good agreement with each other, but neverthe- 
less there are quantitative differences worth remarking 
upon. The baryon density wb becomes much more tightly 
constrained, and looks rather high when compared to 
nucleosynthesis constraints. Perhaps most interestingly, 
assuming HZ the constraints on the reionization epoch 
tighten considerably in favour of early reionization, as in 
abandoning slow-roll we have lost the ability to use tilt 
and tensors to partially mimic the effects of reionization. 



B. A convergence criterion 

As long as the HZ model remains a good fit to the 
data then the issue of introducing the spectral index as 
a parameter, let alone the running of the spectral in- 
dex, needs some consideration. Reconstructing the initial 
power spectrum in bands 25, 2(j indicates that it is ade- 
quately fit by a scale-invariant spectrum, and that there 
should be nothing to be gained by adding any further 
power spectrum parameters. Indeed, why do we include 
further power spectrum parameters at all? If we consider 
the power spectrum as a Taylor expansion 



\nV{k) = \nV{K) + {n s -l)ln( — 



- 1 * 2 (r 

2 V /c* 



(36) 



then the primary reason for including higher power spec- 
trum parameters is to the test the convergence of the 
observable, in this case InV(k), as described at lower- 
order. The first term in the scries has been determined 
to be lnV(k^) ~ —20. We can introduce the next or- 
der parameter, the spectral index, without disrupting the 
constraints on the other parameters, and moreover feel 
confident of our measurement of ln"P(A:*), provided 



\lnV(k*)\ > 



(ng - 1) In 



(37) 



The above criterion is indeed satisfied by all the mod- 
els under consideration since max | ln(fc//c*)| ~ 4 and 
max \n s — 1| ~ 0.07. Thus the determination of the am- 
plitude of scalar perturbations to a reasonable accuracy 
is not questioned by anybody. Similarly, one can feel 
confident in a measurement of n§ if the criterion 



\n s - II > 



(38) 



is satisfied, and the fact that this inequality is typically 
violated for many of the strong running models under 
consideration in the literature {eg. the best fit of Ref. 
is 7is — f = —0.07, as = —0.031 and the two terms are of 
the same order of magnitude) has two possible explana- 
tions. The first is the intriguing possibility that the third 
term in the Taylor expansion Eq. (l.'>(il) dominates over the 
second term at around a scale of ln(fc/fc*) ~ 4, physically 
corresponding to a power spectrum with a maximum near 
fc*. This can only be verified if the contribution to the 
power spectrum from the fourth term (the running of the 
running) is found to be less than the term due to the 
running itself at this scale. The second, and more likely, 
explanation is that the large allowed variation in the run- 
ning of the spectral index is just a symptom of the fact 
that we haven't convincingly determined the spectral in- 
dex yet. We suggest that this criterion can be easily used 
as a check in the fitting procedure: if Eq. I|38|l is violated 
for most models and the running is detected only at low 
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ns 
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BF 


0.023 


0.117 


71.3 


13.1 


23.4 


0.008 


-0.0007 


0.98 


0.07 


BFT 


0.022 


0.115 


70.5 


15.0 


24.1 





0.03 


0.97 





BFHZ 


0.024 


0.119 


71.9 


18.5 


26.8 








1 





\<t> A 


0.022 


0.107 


71.5 


7.1 


20.3 


0.017 


0.017 


0.95 


0.13 


m 2 ct> 2 


0.023 


0.114 


70.9 


10.5 


22.1 


0.008 


0.017 


0.97 


0.06 


TDlo- 


0.025 


0.107 


77.4 


15.0 


23.0 


0.023 


-0.077 


1.03 


0.21 


Bla 


0.025 


0.113 


76.7 


17.2 


24.7 


0.018 


-0.077 


1.04 


0.16 



TABLE I: Cosmological parameters for various inflationary 
models selected by their values of e\ and £2, with power spec- 
trum parameters defined at fc* = O.OlMpc -1 . The models are 
the overall best-fit (BF), the best-fit tilted (BFT), Harrison- 
Zel'dovich (BFHZ), the X(j> 4 and m 2 2 60 e-fold models, the 
best-fit model at the tip of the la contour along the tensor 
degeneracy (TDIct), and the best-fit model with the bluest 
scalar spectrum along the la contour (Blcr). 

significance, then the simplest interpretation is that no 
useful information is coming from the inclusion of the 
new parameter in the fit. One should at least be aware 
of a lack of convergence in the power spectrum observ- 
able, Eq. dnni - 

Alternatively, one can think of Eq. I|38f) as setting the 
boundary between weak and strong running in an obser- 
vational sense. For current observations weak running is 
given by models with 0:3 1 < max|ns — l|/2 ~ 0.03, and 
this can be useful one wishes to introduce a weak running 
prior into the fitting procedure as a perturbation analysis 
to the existing constraints. 

This type of convergence criterion should find echoes 
in other areas of cosmology where an observable is in 
some way series expanded, for instance the dark energy 
equation of state w(z) |27j . 

Before concluding, we bring together the cosmological 
parameters of various inflationary models in Tabled The 
models have been selected by fixing their values of ei and 
£2 , and where a chain was not already available we have 
run a short chain to determine reasonable values for the 
other parameters. These models may be useful as fiducial 
parameters for future study and, being within the la 
contour of Fig. [21 they are approximately degenerate at 
the level of the current data set. 



duced and the underlying inflationary degeneracy now 
becomes visible. Interestingly, if we combine these con- 
straints with a constraint on the number of e-folds of in- 
flation since horizon scale crossing, then we find that the 
A</> 4 inflation model is under strong pressure, though not 
yet definitively ruled out. The m 2 (f) 2 model will come 
under the same threat as long as the data continue to 
favour the Harrison-Zel'dovich spectrum. 

We also introduced a simple convergence criterion, 
Eq. 1)38(1 . to judge the necessity of including higher power 
spectrum parameters such as the spectral index and the 
running of the spectral index. Applying this criterion we 
find that while it is justified to include the spectral in- 
dex in the fit (reflecting the fact that amplitude of scalar 
perturbations is now well determined), it is not useful to 
include the running of the spectral index at present. 

Our inflation analysis comes after those of Refs. [31If3. 
Il3| . with whom we find general agreement. We can make 
the most direct comparison with Barger et al. |12j . who 
used the WMAP data alone with a top-hat Hq prior, and 
a grid based % 2 maximization procedure to obtain their 
constraints. We have have found tighter constraints us- 
ing WMAP, VSA, CBI, ACBAR, 2dF datasets, and an 
MCMC technique, with our 2er constraint in Fig.[2]being 
only slightly looser than their equivalent ler constraint. 
The comparison with Peiris et al. is less straightfor- 
ward, since we did not consider the effect of strong run- 
ning of the spectral index, arguing in Sec. HVBl that we 
can obtain better value from the current data set without 
it. As a result, their results marginally favour a model 
with a blue scalar spectrum on the largest scales (even in 
the tensorless limit), whereas the power-law fits in the lit- 
erature, including their adiabatic/isocurvature fit, favour 
a slightly red spectrum. Both Peiris et al. and Kinney et 
al. [l^l investigate the Monte-Carlo flow reconstruction 
technique, which inevitably involves an extrapolation of 
the inflationary potential and physics well beyond the 
region directly constrained by observations. Therefore a 
certain caution is required when interpreting those anal- 
yses. We have hinted in Sec. IIII CI how one might make 
such an extrapolation of the horizon-flow parameters (in- 
stead of the potential), but the method would require 
considerably tighter constraints on the running of the 
spectral index than we have at present. 



V. CONCLUSIONS 

We have investigated the constraints on various slow- 
roll inflationary models coming from observations of the 
CMB and large-scale structure. The main result is that 
the viable slow-roll parameter space is dramatically re- 
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